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ABSTRACT 

Using new XMM and Chandra observations we present an analysis of the temperature structure of the 
hot gas within a radius of 100 kpc of the bright nearby galaxy group NGC 5044. A spectral deprojection 
analysis of data extracted from circular annuli reveals that a two-temperature model (2T) of the hot 
gas is favored over single-phase or cooling flow (M = 4.5 ± 0.2 M Q yr _1 ) models within the central 
~ 30 kpc. Alternatively, the data can be fit equally well if the temperature within each spherical shell 
varies continuously from ~ Th to T c ~ Tk/2, but no lower. The high spatial resolution of the Chandra 
data allows us to determine that the temperature excursion Th — > T c required in each shell exceeds 
the temperature range between the boundaries of the same shell in the best-fitting single-phase model. 
This is strong evidence for a multiphase gas having a limited temperature range. We do not find any 
evidence that azimuthal temperature variations within each annulus on the sky can account for the range 
in temperatures within each shell. We provide a detailed investigation of the systematic errors on the 
derived spectral models considering the effects of calibration, plasma codes, bandwidth, variable Ah, and 
background rate. We find that the RGS gratings and the EPIC and ACIS CCDs give fully consistent 
results when the same models are fitted over the same energy ranges for each instrument. The cooler 
J> ■ component of the 2T model has a temperature (T c ~ 0.7 keV) similar to the kinetic temperature of the 

' stars. The hot phase has a temperature (Th ~ 1.4 keV) characteristic of the virial temperature of the 

, ~ 10 13 Mq halo expected in the NGC 5044 group. However, in view of the morphological disturbances 

and X-ray holes visible in the Chandra image within R « 10 kpc, bubbles of gas heated to ~ Th in this 
region may be formed by intermittent AGN feedback. Some additional heating at larger radii may be 
associated with the evolution of the cold front near R ~ 50 kpc, as suggested by the sharp edge in the 
EPIC images. 

Subject headings: X-rays: galaxies: clusters - galaxies: halos - galaxies: formation - cooling flows - 
galaxies: individual: NGC 5044 
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r . r .. i 4.1 and spectral properties of the diffuse hot gas out to radii 
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,,,,,, ^, • u n hot gas on smaller scales. 
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, t ... (1 , . , j, ■ In tins paper we analyze tnc temperature structure ot 

suited for studies ot their temperature structure since their ,. . , ^ " , T „„ .„;, _, , , , . , 

, , it . 4., , . ... the hot gas m NGC 5044. Ihe metal abundances Buote 

~ 1 kev temperatures occur near the peak m sensitiv- . , r . n i ■. .■ i- 

., c . v j 4-1, 4- i • ii et al. 2003b, hereafter Paper 2) and gravitating mass dis- 

lty ot current X-ray CCDs and near the strong, highly ■,• ,• , . • 

... ■ t u ii • • t tribution are discussed m companion papers, 
temperature-sensitive iron L-snell emission lines. 1 1 r 

Using new XMM and Chandra observations we present 2 observations and data preparation 

an analysis of the temperature structure of the hot gas 

of NGC 5044, perhaps the brightest group in soft X- 2.1. XMM 

rays. Although several previous ROSAT and ASCA stud- NGC 5044 was observed with the EPIC pn and MOS 

ies demonstrated that the hot gas within r - 30 kpc is not CCD cameras for approximately 20 ks and 22 ks respec- 

isothermal (e.g., David et al. 1994; Buote 1999), the data t i ve i y during AO-1 as part of the XMM Guest Observer 

could not distinguish between single-phase and multiphase program. We generated calibrated events lists for the data 

models. The combined spatial and spectral resolution of using tne standard SAS v5.3.3 software. Since the diffuse 

the XMM and Chandra CCDs allows for unprecedented emission of NGC 5044 fills the entire field of view we es- 

mapping of the temperatures and elemental abundances timate the background using the standard "background 

of the hot gas in groups and clusters. The higher energy templates". These templates are events lists obtained by 
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combining several high Galactic latitude pointings 4 . 

Inspection of the CCD light curves for events with ener- 
gies above 10 keV does not reveal any strong flares in the 
NGC 5044 observation, but there are times of increased 
activity which are more pronounced in the pn data. After 
applying the count-rate screening criteria recommended to 
match the background templates (and also the standard 
screening criteria recommended for all observations), we 
arrive at final exposures of 19.5 ks for the MOS1, 19.3 ks 
for the MOS2, and 8.9 ks for the pn. 

Since the quiescent background varies typically by ~ 
10% it is necessary to normalize the background tem- 
plates to each source observation. Since NGC 5044 has 
a gas temperature ~ 1 keV there is little emission from 
hot gas for energies > 5 keV. We renormalized the MOS1 
and MOS2 background templates by comparing source and 
background counts in the 7-12 keV band extracted from 
regions near the edges of the MOS fields. We obtain back- 
ground normalizations that are %16 and %5 above nominal 
respectively for the MOS1 and MOS2. 

Since the strict events screening mentioned above leaves 
the pn with < 50% of its raw exposure, it is evident that 
background flares contaminate the pn more than the MOS. 
Even after the strict screening, we still require a back- 
ground normalization that is 21% above nominal for the 
pn. Although this fraction is not much larger than for 
the MOS1, an excess of the source over the background is 
clearly visible above energies of a few keV. We have there- 
fore also investigated a potentially more accurate method 
of subtracting the flaring background following our study 
of the XMM observation of NGC 1399 (Buote et al. 2003a). 
That is, after subtracting the pn spectra taken from re- 
gions near the edge of the field with the standard back- 
ground templates re-scaled by their nominal exposures, 
we fitted the resultant spectrum with a two-component 
model consisting of a thermal component, represented by 
an APEC thermal plasma, and a broken power-law (BPL) 
model, representing the residual flaring background which 
is most pronounced at high energies. For the BPL model 
we obtain a break energy of 0.6 keV with power-law in- 
dices of -2.6 below the break and 0.54 above the break; 
i.e., over most of the bandpass of interest (> 0.6 keV) the 
flaring background is well-described by a power-law with 
index 0.54. 

This BPL model defines the shape of the excess back- 
ground above the standard templates. We determined its 
normalization separately for each region of interest using 
data between 6.0 and 7.25 keV. This 6 keV lower limit is 
selected to avoid contamination from softer source emis- 
sion while the upper limit is chosen to avoid calibration 
emission lines. We find that this method to subtract the 
excess background in the pn data gives results consistent 
with the simpler method of renormalizing the background 
template. In this paper we use the BPL model of the ex- 
cess background as our default method for the pn. 

2.2. Chandra 

4 See XMM calibration note by D. Lumb (XMM-SOC-CAL-TN-0016). 

5 http:/ /cxc. harvard.edu/udocs/docs/docs. html 

6 http://cxc.harvard.edu/ciao/ 

7 http://cxc.harvard.edu/cal 



NGC 5044 was observed by Chandra with the ACIS-S3 
camera for w 22 ks during AO-1. The events list was cor- 
rected for charge-transfer inefficiency according to Towns- 
ley et al. (2002), and only events characterized by the stan- 
dard ASCA grades 5 were used. The standard CIAO 6 soft- 
ware (version 2.2.1) was used for most of the subsequent 
data preparation. 

Since the diffuse X-ray emission of NGC 5044 fills the en- 
tire S3 chip, we used the standard background templates 7 
to model the background. After running the standard 
LC_CLEAN script to clean the source events list of flares 
with the same screening criteria as the background tem- 
plates, we arrive at a final exposure time of 20.2 ks for 
NGC 5044. 

3. IMAGE AND RADIAL PROFILE 

In Figure 1 we display an adaptively smoothed mosaic 
of the XMM MOS1 and MOS2 images as well as the adap- 
tively smoothed Chandra ACIS-S3 image. Several sources 
can be seen embedded in the diffuse emission which fills 
the entire MOS field. That we see more sources at large 
radius suggests that discrete sources at smaller radii are 
swamped by the diffuse emission which is brightest near 
the center. A careful inspection of the MOS and pn im- 
ages within R as 1' of NGC 5044 reveals significant non- 
circular features. The higher-resolution Chandra ACIS-S 
image (Figure 1) shows these features to be holes in the hot 
gas similar to those seen in Chandra images of the cores 
of galaxy clusters (e.g., Fabian et al. 2000; McNamara 
et al. 2000; Churazov et al. 2001). Typically these holes 
are filled by extended radio emission that emanates from 
a central source. However, the (a 45") NVSS (Condon 
et al. 1998) radio image of NGC 5044 reveals only a fairly 
bright central point-like source suggestive of an AGN. Per- 
haps a deeper observation at higher resolution would find 
extended radio emission filling the X-ray holes. 

Over a radius R w 1 — 5' (i.e., to the edge of the central 
MOS CCD) the diffuse emission is approximately circu- 
lar with concentric isophotes. Between R w 5 — 6' the 
isophote centroid is slightly offset from values at smaller 
radius. The origin of this shift is apparent in Figure 1 
as excess emission to the East and N-E outside the cen- 
tral CCD. This offset was noted previously in the ROSAT 
PSPC observation of NGC 5044 (David et al. 1994). The 
XMM image, however, also reveals that there is a sharp 
edge in the X-ray surface brightness on the West, N-W 
side near this offset region, R <~ 6'. This type of sharp 
feature has been observed recently in Chandra observa- 
tions of several galaxy clusters and has been interpreted 
as a "cold front"; i.e., a contact discontinuity between re- 
gions of different density and temperature (see Markevitch 
et al. 2002). We note that at radii larger than ~ 6' the 
isophote centroids are consistent with those for R<5'. 

We characterized the image further by computing the 
azimuthally averaged surface brightness profile in the 0.3- 
5 keV band separately for the MOS and pn data. We 
created exposure maps for each detector using the stan- 
dard SAS software and used these maps to flatten the 
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Fig. 1. — (Left panel) False-color mosaic of the XMM MOS1 and MOS2 images of NGC 5044 adaptively smoothed using the SAS task 
ASMOOTH. The image has been divided by the summed MOS1 and MOS2 exposure maps to correct for exposure variations. (Right panel) 
Portion of the archival Chandra ACIS-S image (also false-color) located at the center of NGC 5044 adaptively smoothed using the CIAO task 
CSMOOTH. In both images celestial N is up and E is to the left. 
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Fig. 2. — (Left panel) Azimuthally averaged 0.3-5 keV surface brightness profile of the background-subtracted, exposure-corrected combined 
MOS1 and MOS2 data (kpc units on top axis). The lower points represent the combined background values. The best-fitting /3 model for the 
MOS data is indicated by the dashed line. (Right panel) Same quantities are plotted for the pn data. 



images. We combined the MOS1 and MOS2 data into 
a single exposure-corrected image. The same procedures 
were followed for the background images. (Since here we 
are not concerned with the detailed spectral dependence of 
the background for the pn we used the renormalized tem- 
plate accounting for the factor of 21% discussed in §2.) 
Point sources and chip gaps were masked out from the 
calculation. 



The radial profiles, £ x (r), of the source and background 
images are shown in Figure 2. Within R ~ 7' the radial 
profiles are similar for both detectors, but at larger radii 
the profiles appear to deviate. To facilitate comparison 
between the profiles, we fitted standard "/3 models" con- 
volved with the XMM PSF 8 to S x (r). The best-fitting 
models are shown in Figure 2. We obtained parameters, 



See XMM calibration note EPIC-MCT-TN-012 by Ghizzardi. 
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r c = 61.9 ± 0.4", and (3 = 0.585 ± 0.001 for the MOS and, 
r c = 47.9 ± 0.4", and (3 = 0.522 ± 0.001 for the pn. The 
values of these parameters confirm the impression from the 
figure that at large radius S x (r) is steeper for the MOS. 

The steeper profile of the MOS at larger radii is proba- 
bly the result of small errors in the background estimate at 
large radius. Near R = 7' the MOS background rate equals 
the source rate indicating that the derived S x (r) for the 
source is very sensitive to small errors in the background 
for R > 7'. In contrast, the background rate does not 
equal the source rate for the pn until near the edge of the 
field, R <~ 11'. The robustness of £ x (f) for the pn is also 
supported by the analysis of the ROSAT data of NGC 5044 
by David et al. (1995) who obtained (3 = 0.53 ± 0.02 con- 
sistent with our result for the pn. If data with R > 400" 
are excluded for the MOS, then we obtain best-fitting pa- 
rameters, r c = 52.1" and j3 = 0.542, in better agreement 
with the pn and ROSAT values. 

We conclude that a single (3 model is a good average 
description of the radial surface brightness profile of NGC 
5044 between R ~ 10" and ~ 11'. (At the center there 
is evidence for a small excess above the (3 model on the 
scale of the PSF. This is also the region where the az- 
imuthal distortions are most pronounced in the Chandra, 
image.) Therefore, the possible "cold front" suggested by 
the isophotal centroid offset near R = 6' does not sub- 
stantially distort the radial profile away from that of an 
equilibrium configuration. 

4. SPECTRAL DEPROJECTION ANALYSIS 

Motivated by the regular appearance of the radial pro- 
file and the lack of substantial azimuthal variations in 
the spectral properties, we focus our analysis on the az- 
imuthally (and spherically) averaged spectral properties 
of the X-ray emission. We present an analysis of the az- 
imuthal spectral properties of the Chandra and XMM data 
in §5. 

4.1. Preliminaries 

For each XMM detector (MOSl,MOS2,pn) we extracted 
spectra in concentric circular annuli located at the X-ray 
centroid of the inner contours (i.e., computed within a 
2' radius) such that the width of each annulus contained 
~ 8000 background-subtracted counts in the 0.3-5 keV 
band in each MOS detector, and the minimum width was 
set to 1' for PSF considerations; note that we address 
residual effects of the PSF by comparison to fits of the 
high-resolution Chandra data alone in §4.3 and §6. Ob- 
vious point sources were masked out before the extrac- 
tion (and corresponding regions also were omitted from 
the background templates). These restrictions resulted in 
a total of eight annuli within r = 10'; these annuli are de- 
fined in Table 1. (Note that the results we present below 
are quite insensitive to the choice of radius within which 
to compute the centroid.) RMF and ARF files were gen- 
erated using SAS. (Note that the vignetting correction is 
administered through the ARF file.) Finally, the spectral 
pulse-invariant (PI) files for each annulus were rebinncd 
such that each energy bin contains a minimum of 30 counts 
appropriate for \ 2 fitting. 

We first consider the Chandra data extracted in the 
same regions as the XMM data subject to the constraint 



that the regions lie entirely within the ACIS-S3 field. This 
restriction gives Chandra data within the first three radial 
bins defined for the XMM data. Later in §4.3 we shall ex- 
amine the Chandra data within smaller spatial regions. 
We have applied the latest corrections to the Chandra 
ARF files that account for a time-dependent degradation 
in the quantum efficiency at lower energies (i.e., CIAO COR- 
rarf routine). 

To obtain the three-dimensional properties of the X-ray 
emitting gas we perform a spectral deprojection analysis 
assuming spherical symmetry using the (non-parametric) 
"onion-peeling" technique. That is, one begins by deter- 
mining the spectral model in the bounding annulus and 
then works inward by subtracting off the spectral contri- 
butions from the outer annuli. We have previously devel- 
oped a code based on xspec (Arnaud 1996) to implement 
this procedure (Buote 2000a). Consequently, we obtain 
temperatures, abundances, densities, and absorption col- 
umn densities (and any other desired parameters) by fit- 
ting spectral models to the "deprojected spectra" . 

In this procedure we account for the emission projected 
from shells outside of our bounding annulus by assuming 
the X-ray emissivity for such shells varies as a power-law 
with radius and has the same spectral shape as deter- 
mined from the bounding annulus. For the XMM data 
for NGC 5044 we take the radial surface brightness to 
vary as r -2 (or /3 = 0.5 for the /3 model) as indicated 
by the data within R = 10'. As found previously (§3.1 
of Buote 2000a) the derived spectral parameters (except 
density in the penultimate annulus) are quite insensitive 
to this choice. 

To estimate the uncertainties on the fitted parameters 
we simulated spectra for each annulus using the best- 
fitting models and fit the simulated spectra in exactly the 
same manner as done for the actual data. From 20 Monte 
Carlo simulations we compute the standard deviation for 
each free parameter which we quote as the "lcr" error. 
(This is a slightly different manner of quoting the errors 
compared to §3.2 of Buote 2000a.) 

Because deprojection always inflates the errors between 
points, which is related to the error associated with the 
derivative of the emissivity in an Abel inversion, we reg- 
ularize (i.e., smooth) some of the parameters (see §3.3 of 
Buote 2000a) . In actuality we regularize ex post facto by 
restricting the value of a parameter within certain bounds 
specified by a pre-determined radial variation for the pa- 
rameter. We find it necessary to regularize only the O and 
Ne abundances so that the radial abundance variation of 
the logarithmic derivative is between ±2. We emphasize 
that rcgularization only applies to O and Ne for the depro- 
jected (3D) spectral analysis. No regularization is applied 
to any 2D model. 

We take the solar abundances in xspec (vl 1.2. Oaf) to be 
those given by the Grevesse & Sauval (1998) table which 
use the correct new photospheric value for iron which 
agrees also with the value obtained from solar-system me- 
teorites (e.g., McWilliam 1997). 

4.2. Simultaneous Fitting of XMM- Chandra Data 

Within the eight annuli defined in Table 1 we consider 
simultaneous fits to the MOS1, MOS2, and pn XMM data 
within radial bins 4-8 for which we do not have Chan- 
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Table 1 

Quality of Spectral Fits (x 2 /dof) for IT and 2T Models 







Rout 




IT 


2T 


Bin 


(arcmin) 


(arcmin) 


2D 


3D 


2D 


3D 


1 


0.0 


0.5 


779.8/451 


737.9/451 


646.7/449 


644.0/449 


2 


0.5 


1.5 


1812.3/689 


1368.1/689 


974.2/687 


949.8/687 


3 


1.5 


2.5 


1338.2/693 


1189.0/693 


855.3/691 


836.0/691 


4 


2.5 


3.5 


730.6/482 


666.8/482 


610.3/480 


537.3/480 


5 


3.5 


4.5 


535.6/454 


504.7/454 


521.7/452 


484.4/452 


6 


4.5 


5.7 


547.7/469 


548.3/469 


516.2/467 


514.4/467 


7 


5.7 


7.6 


538.0/515 


572.4/515 


519.7/513 


550.0/513 


8 


7.6 


10.1 


674.6/561 


674.6/561 


658.7/559 


658.7/559 



Note. — Bin refers to an annulus for 2D models and to a shell for 3D models. Single-temperature (IT) models are described in §4.2.1 and 
two-temperature (2T) models are described in §4.2.2. 



dra data. Within radial bins 1-3 wc perform simultaneous 
fits to the Chandra and XMM data. Comparison of the 
Chandra and XMM data within their overlap regions is 
discussed in §6. We focus on fits performed over the en- 
ergy range 0.5-5 keV from consideration of the iron abun- 
dances obtained from separate Chandra and XMM fits as 
discussed in Paper 2. However, all results presented for 
the temperatures in this paper are consistent with those 
obtained when fitting 0.3-5 keV (§6). 

Our baseline model (IT) consists of a single thermal 
plasma component using the APEC code modified by fore- 
ground Galactic absorption (TVh = 5 x 10 20 cm~ 2 ) using 
the PHABS model in xspec. The free parameters for this 
baseline model are all associated with the plasma compo- 
nent: temperature (T), normalization, and Fe, O, Ne, Mg, 
Si, and S abundances - all other elements tied to Fe in 
their solar ratios. (Models with intrinsic absorption are 
discussed in §6.) We obtained results for the IT model 
fitted directly to the data projected on the sky (i.e., tra- 
ditional 2D model) and also fitted to the deprojected data 
(i.e., 3D model) as discussed in §4. 

To account for small calibration differences between de- 
tectors we multiplied the spectrum of each detector by a 
detector-dependent constant. These constants were deter- 
mined by fitting the IT model to the data with the con- 
stants as free parameters. We fixed the constants to these 
values for all subsequent spectral fits. This approach guar- 
antees that the relative normalizations of multiple emission 
components (e.g., in 2T models) are the same for each de- 
tector. 

4.2.1. Single- Temperature Models 

The temperature profiles for the IT models (2D and 
3D) arc displayed in Figure 3. For the IT (2D) model 
the temperature rises from T <~ 0.7 keV at the cen- 
ter to T <~ 1.2 keV at large radii consistent with previ- 
ous 2D ROSAT determinations (e.g., David et al. 1994; 
Buotc 2000a). The 3D temperature has a steeper rise out- 
ward from the center consistent with the 2D profile being 
smeared by projection effects. The lower temperature in 
the bounding annulus certainly suggests a declining tem- 
perature value, but its value is underestimated because of 
projection; i.e., we have had to assume the projected emis- 
sion from shells exterior to our bounding annulus has the 
same spectral shape as the bounding annulus (§4). 



The x 2 values for both models are listed in Table 1. 
This model is a reasonably good fit in the outer annuli, 
but the fit degrades at progressively smaller radii, though 
the trend reverses somewhat for the central annulus. For 
most radial bins the 2D and 3D fits are similar, with the 
3D model usually giving slightly better x 2 values. In shells 
2-4, however, the 3D model clearly provides a better fit 
indicating that annuli 2-4 contain multiple temperature 
components projected on to the sky. 

The IT (3D) model is still a relatively poor fit in shells 
1-4. In Figure 5 we show the EPIC and ACIS spectra for 
annulus 2 fitted to the IT (3D) model fits and residuals. 
The IT (3D) model yields fit residuals near 1 keV that are 
fully characteristic of those obtained when trying to force 
a single-temperature model to fit a spectrum consisting of 
a multiple components with temperatures near 1 keV (e.g., 
Buote & Fabian 1998; Buote 1999, 2000b). The small de- 
viations above 2 keV also suggest the presence of another 
(higher) temperature component. 

4.2.2. Two-Temperature Models 

Since the fits within the central regions - even when the 
spectral data are deprojected - indicate the presence of 
at least one more spectral component, we investigated fit- 
ting simple two-temperature models (2T) to the data. The 
abundances for each component are tied together to limit 
the number of free parameters and because relaxing this 
condition does not amount to a substantial improvement 
in the fit quality. (In Paper 2 we briefly discuss 2T models 
with separately varying abundances on each component.) 
Consequently, the 2T models add only two free param- 
eters: the temperature and normalization of the second 
component. 

The x 2 values for both the 2T (2D) and 2T (3D) mod- 
els are listed in Table 1. In all spatial bins the fits are 
improved by the 2T models, with the greatest improve- 
ment occurring in the central bins (1-4). The considerable 
improvement for bin 2 is shown in Figure 6 for the 2T 
(3D) model. In comparison to Figure 5 the fit residuals 
are reduced substantially, and overall the fit is as good as 
could be expected for such a simple model. 

The values of X 2 obtained for the 2T (2D) and 2T (3D) 
models are very similar, but the 3D models generally have 
lower x 2 values. Only in bin 4 is the 3D model an obvi- 
ous improvement over the 2D version. Since 2T models 
can mimic radial temperature variations accurately, this 
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Fig. 3. — Radial temperature profiles (units - bottom: arcminutes, top: kpc) and la errors for (Left panel) IT and (Right panel) 2T models 
fitted simultaneously to XMM and Chandra data. Note that the Chandra data apply only to the inner three radial bins. In each case "3D" 
refers to results obtained from a spectral deprojection analysis. 
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Fig. 4. — Emission measure of the cooler and hotter components for the 2T (3D) model whose temperatures are plotted in Figure 3. The emis- 
sion measure is expressed as the normalization of the APEC plasma model as given in XSPEC; i.e., normalization is 10 — 14 n e n p V/ (AttD 2 a (1 + z) 2 ), 
where n e is the electron number density, n p is the proton number density, V is the emitting volume, Da is the angular diameter distance, 
and z is the redshift. The units are arcminutes on the bottom axis and kpc on the top axis. 



agreement is not unexpected. Nevertheless, because of the 
similarity between the 2D and 3D fits, projection effects 
are not very important for the 2T models. 

The temperature profiles of the 2D and 3D 2T models 
are displayed in Figure 3. Focusing for the moment on the 
2D model which has better constrained values than the 
3D model, it is seen that both the cool (T c ) and hot (Th) 
temperature components are consistent with an isothermal 
radial profile for R > 20 kpc. For smaller radii, T c remains 
nearly constant with an indication of a slight downturn in 
the central two bins. The hotter component rises slightly 
at the center. The values of T c at small radii and the val- 
ues of Th at large radii are similar to the IT temperature 
values. 

Within the 1-2 a errors the 3D and 2D temperatures 



are consistent in all bins, but there are some trends that 
deserve attention. At large radii (r > 40 kpc) T c is not 
tightly constrained particularly for the 3D model; the ex- 
ception is the final bin where the 2D and 3D fits are equiva- 
lent in our deprojection procedure (apart from the inferred 
normalization). In bin 7 the cooler component is barely de- 
tected in 3D; the ratio of normalizations for the best-fitting 
cooler-to-hotter components is 0.12 for the 2D model com- 
pared to only 0.03 for the 3D model. Also in bin 7 the value 
of T c is poorly constrained for the 3D model, and the best- 
fitting value of 0.5 keV is located at the lower range we 
allowed for during the fit. We attribute the stronger de- 
tection of the 2T model in bin 8 compared to bin 7 to 
the limitation in the way we treat the deprojection of the 
bounding annulus; i.e., our deprojection method assumes 
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the source flux at large radii has the same spectrum as the 
bounding shell. 

At small radius (r < 10 kpc), the temperature of the 
hotter component rises and becomes poorly constrained 
in bin 1 for the 3D model, where the best-fitting value 
isT = 3.7 keV and the lowest value from 100 Monte Carlo 
simulations is 1.4 keV which is shown in Figure 3. It is 
possible that the rise in Th in bin 2, and especially bin 
1, arises from the emission of unresolved discrete sources 
similar to the trend observed for the related system, NGC 
1399 (Buote 2002). If we add a 10 keV bremsstrahlung 
(brcm) component to the 2T models the fits are not im- 
proved, but the fitted value of Th in bins 1-2 are consistent 
with those adjacent bins. Moreover, we obtain a luminos- 
ity L x « 3 x 10 40 erg cm~ 2 s _1 for the brcm component. 
Since this value is comparable to the (relatively uncertain) 
luminosity expected from discrete sources in NGC 5044 
using the results of O'Sullivan et al. (2001), we conclude 
that the unresolved emission from discrete sources is a vi- 
able candidate for the rise in Th at small radius in the 2D 
models (though Th is still a separate component from the 
discrete sources). 

In Figure 4 we plot the emission measures for each com- 
ponent of the 2T (3D) model. The cooler component dom- 
inates within ~ 10 kpc while the hotter component dom- 
inates for r > 50 kpc. Over w 15 — 40 kpc the relative 
emission measures of the components are within a factor 
of two of each other. 



to this cooling flow model as "CF+1T". Note that this 
model adds only one free parameter over the IT case, i.e., 
M, since we tie together the abundances between the two 
components. 

The x 2 values for the CF+1T model fitted to the XMM 
and Chandra data are listed in Table 2. Although the 
CF+1T fits are an improvement over the IT models in 
some bins, the improvement is not nearly as great as ob- 
served for the 2T model. The largest improvement occurs 
in shells 3-4 where the \ 2 values of the CF+1T model 
are roughly half-way between the IT and 2T values. The 
CF+1T model offers no improvement over the IT model 
in bins 1-2. 

The total mass deposition rate (see Table 3) for all shells 
(i.e., within r = 97 kpc) is, M = 4.5 ± 0.2 M Q yr _1 . 
This value is less than half that obtained from the ROSAT 
imaging analysis of NGC 5044 by David et al. (1994) and 
only « 1/8 that obtained from the single-aperture analy- 
sis of the ASCA data by Buote (1999). The much smaller 
value of M obtained from the spatially resolved spectral 
analysis of the XMM and Chandra data is very similar to 
the reduction obtained from Chandra and XMM results 
for cooling flow clusters (e.g., David et al. 2001). 

In fact, M is even slightly smaller because the value in 
the bounding shell is very likely overestimated. The spike 
in M in shell 8 is almost certainly biased by the edge effect 
inherent in our deprojection procedure (§4). 



4.2.3. Cooling Flow and Other Multitemperature Models 

Since the 2T models are preferred over the IT models 
(particularly in bins 1-4) , we investigated other multitem- 
perature spectral models to attempt to constrain the gen- 
eral differential emission measure (DEM) of the emission 
spectrum. Because the task of inferring the general DEM 
from the X-ray spectrum of a coronal plasma is difficult 
and often highly degenerate (Craig & Brown 1976), espe- 
cially for data of only moderate energy resolution such as 
provided by the XMM and Chandra CCDs, we probed the 
DEM using a few parameterized models. For this discus- 
sion we focus on models fitted to the deprojected spectra; 
i.e., 3D models. 

Cooling Flow: The DEM for an ideal gas cooling at con- 
stant pressure from a maximum temperature, T max , is 
(e.g., Johnstone et al. 1992), 

UT) = ^ 



2 fim p A(T) ' ^ 

for T < T max . Here ks is Boltzmann's constant, m p is the 
proton mass, A(T) is the plasma emissivity which we take 
to be given by the APEC code, and M is the mass drop-out 
rate. The luminosity of this cooling gas emitted over some 
energy range AT is, 



/■Tr 
cf _ / 

JAE ~ T ■ 

u -L m r 



A AE (T)^ ci (T)dT, 



(2) 



where Aas(T) is the plasma emissivity integrated over the 
bandwidth AT. We set T min = 0.05 keV. The total emis- 
sion spectrum of a cooling flow model is taken to be L C ^ E 
plus the spectrum of a IT model with T = T max to rep- 
resent the emission from the ambient gas. We shall refer 



Gaussian: An important consideration is whether the 
better fits obtained by the 2T (3D) model over the IT 
(3D) model could simply reflect the range of temperatures 
within each shell arising from the radially varying temper- 
ature profile implied by the IT fits. In support of this hy- 
pothesis is the fact that the improvement of the 2T model 
is greatest in shells 1-4 where the IT temperature profile 
is changing most rapidly with radius; i.e., there we should 
have the largest range of temperatures per shell. 

To test this hypothesis we have explored a model where 
the temperature distributions within each shell are a Gaus- 
sian. The Gaussian DEM (GDEM) is expressed as, 



4 CO 



(TTV 2-7T 



exp [-(T-T ) 2 /24] , 



(3) 



where T is the mean and ot the standard deviation of the 
Gaussian, and norm is a constant proportional to n e n p V 
as defined for the APEC code in xspec (see caption to 
Figure 4). The luminosity over some energy range AT is 
therefore given by equation (2) with £ g (T) replacing £ c f (T) 
and T max and T m i n set to ±oo. 

In shells 1-4 the GDEM model provides superior fits to 
the IT and cooling flow models, but does not have x 2 va l~ 
ues quite as low as the 2T model. For example, in shell 3 
we obtain for the Gaussian model, x 2 = 869.7 for 692 dof 
compared to X 2 = U89.0 for IT, X 2 = 982.1 for CF+1T, 
and x 2 — 836.0 for 2T. Note that the Gaussian adds only 
ut as a free parameter over the IT model; i.e., it has the 
same number of free parameters as the cooling flow model 
but one less than the 2T model. 

The fitted values of o~t for shells 5-7 are consistent with 
the single-phase hypothesis within the 1 — 1.5(7 errors on 
ct t . For shell 8 a T = 0.22 ± 0.05 keV could reflect the 
projection of cooler gas from exterior shells as discussed 
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Table 2 

Quality of Spectral Fits (x 2 /dof) for other Multitemperature Models 



Shell 


Cooling Flow 


GDEM 


PLDEM 


1 


729.6/450 


695.5/450 


665.0/449 


2 


1375.6/688 


1091.4/688 


973.3/687 


3 


982.1/692 


869.7/692 


847.2/691 


4 


591.2/481 


546.5/481 


546.9/480 


5 


490.7/453 


484.6/453 


494.7/452 


6 


526.0/468 


518.9/468 


535.8/467 


7 


561.1/514 


554.3/514 


558.1/514 


8 


659.7/560 


658.3/560 


660.6/560 



Note. — All models arc 3D. See §4.2.3 for description of the models. 



in §4.2.2 for the 2T models. We cannot determine this for 
certain with the present data set. 

In contrast, in shells 2-4 (and perhaps shell 1) the fit- 
ted values of To and ot do not appear to be consistent 
with the radially varying single-phase hypothesis. The lor 
ranges for the temperatures of the GDEM model arc 0.68- 
0.97 keV for shell 2, 0.78-1.11 keV for shell 3, 0.95-1.48 keV 
for shell 4. The considerable overlap of these lcr values 
is in conflict with the small temperature range expected 
from the single-phase hypothesis. That is, for the IT (3D) 



model, the temperature difference between shells 1 and 3 
is 0.21±0.01 keV, between shells 2 and 4 is 0.32±0.02 keV, 
and between shells 3 and 5 is 0.35 ± 0.02 keV. The single- 
phase hypothesis predicts, therefore, a range of tempera- 
tures of approximately half the difference between adjacent 
shells; i.e., ss 0.11 keV within shell 2, ss 0.16 keV within 
shell 3, and « 0.17 keV within shell 4. These values are 
considerably smaller than the \ot ranges of 0.28±0.02 keV 
for shell 2, 0.32 ± 0.02 keV for shell 3, and 0.52 ± 0.06 keV 
for shell 4. 



9 



Table 3 

Selected Parameters for other Multitemperature Models 





(!r»(~il i vi cr Flow 

WUUllllg J_ 1UW 


GDEM 




PLDEM 






M 


T 


(7 T 


a 


T ■ 

■*■ min 


T — T ■ 

± max ± min 


Shell 


(M e yr" 1 ) 


(kcV) 


(keV) 




(kcV) 


(keV) 


1 


0.00 ±0.15 


0.72 ± 0.01 


0.10 ±0.04 


-5.0 ±2.3 


0.62 ±0.14 


> 0.5 


2 


3.02 ±0.13 


0.83 ±0.01 


0.14 ±0.01 


-4.6 ±0.5 


0.69 ±0.01 


> 0.7 


3 


0.71 ±0.07 


0.94 ±0.01 


0.16 ±0.01 


-3.9 ±2.1 


0.75 ± 0.04 


1.2 ± 0.6 


4 


0.33 ±0.05 


1.22 ±0.02 


0.26 ± 0.03 


0.6 ± 1.6 


0.72 ±0.13 


1.0 ± 0.2 


5 


0.01 ±0.03 


1.31 ±0.05 


0.19 ±0.11 


3.4 ±4.8 


1.19 ±0.30 


0.1 ± 0.6 


6 


0.14 ±0.04 


1.37 ±0.06 


0.29 ±0.14 


-5.4 ±3.2 


1.19 ±0.28 


0.1 ± 0.6 


7 


0.04 ±0.07 


1.27 ±0.03 


0.00 ± 0.04 





1.23 ±0.35 


0.1 ± 0.8 


8 


0.22 ±0.05 


1.22 ± 0.02 


0.22 ±0.04 





0.87 ±0.27 


0.7 ±0.7 



Note. — All models arc 3D. See §4.2.3 for details of the models. When a lower limit is given it represents the lowest value obtained from 20 
error simulations; i.e., it is essentially a 95% confidence lower limit. 



The temperature difference between shells f-2 is 0.09 ± 
O.Of keV for the IT (3D) model implying a single-phase 
temperature range of ~ 0.05 keV. The Ictt range for the 
GDEM model in shell 1 is a T = 0.10 ± 0.04 keV which is 
only marginally discrepant with the single-phase hypoth- 
esis. 

Power Law: Finally, to further quantify deviations from 
the single-phase hypothesis we have investigated spectral 
fits using a power-law DEM (PLDEM), 

{norm rp a -i 

log(T max /T mm ) J , a--L 

norm( a +l) rra n / , W 
(T°+ X 1 -T Q + 1 ) ' ^ 

where norm is a constant proportional to n e n p V as de- 
fined for the APEC code in XSPECas discussed previously 
(see caption to Figure 4. The luminosity over some energy 
range AE is therefore given by equation (2) with £ p i(T) 
replacing £ c f(T). This model adds two free parameters 
over the IT model: a and the width of the temperature 
distribution, T max — T m i n ; i.e., the power-law model has 
the same number of free parameters as the 2T model. Be- 
cause the values of a and T max — T min are correlated we 
found it necessary to restrict their values during the fits 
to -5 < a < 5 and T max - T min < 10 kcV. 

As indicated in Tabic 2, the PLDEM model fits nearly as 
well as the 2T model in every shell. Because a is not well 
constrained, particularly at larger radii, we fixed a = 
in in shells 7-8. Even with this restriction, the power-law 
model fits nearly as well as the 2T model in those shells. 

The results obtained for a, T m i n , and T max — T m i n are 
listed in Table 3. Like both the 2T and GDEM models, 
in shells 5-7 there is little indication of significant multi- 
temperature gas; i.e., T m - ln > 1.2 keV and T max — T min is 
small (though uncertain). Similarly, in shell 8 the PLDEM 
model indicates the presence of multiple temperature com- 
ponents which likely arises from projection of gas exterior 
to shell 8 (see §4.2.2). 

For shells 1-4 the PLDEM fits imply a substantial range 
of temperature components in each shell. In shell 4 we 
have a w which indicates an equal contribution from 
temperature components over a range from approximately 
0.7 keV to 1.7 keV (with some uncertainty in the upper 
limit). The value of a becomes increasingly negative at 



smaller radius showing that the temperature distribution 
is progressively more peaked near r m ; n . However, the large 
values of T max — T m i n reflect the rise in Th seen for the 2T 
models at lower radius which is probably attributed to the 
emission from unresolved stellar sources (§4.2.2). 

Overall, the PLDEM model, like the 2T and GDEM 
models, indicate that within shells 1-4 the temperature 
distribution is too wide within each shell to be described 
by a radially varying single-phase medium. The XMM 
CCD data cannot distinguish between the PLDEM and 
2T models; i.e., a continuous versus a discrete tempera- 
ture distribution. 

4.3. Smaller Regions with Chandra 

The higher spatial resolution of Chandra allows for 
smaller regions to be probed than with XMM. This allows 
the preference for multitemperature models over the ra- 
dially varying single-phase hypothesis of NGC 5044 to be 
tested with Chandra using shells of smaller width. If the 
hot gas is really a radially varying single-phase medium, 
then the temperature range indicated by a 2T model 
should decrease with decreasing aperture size. Therefore, 
we have further divided the region spanned by shells 1- 
3 into the six shells 1C-6C defined in Table 4. We are 
only able to usefully probe apertures of about half the size 
used previously because the Chandra ACIS-S3 data alone 
do not have the combined sensitivity of the XMM CCDs. 

In Tabic 4 we list the \ 2 values for the IT and 2T (both 
3D) fits in shells 1C-6C. In good agreement with the re- 
sults presented above, we find that the 2T model is a better 
fit than the IT model. The temperatures of the IT and 
2T models are plotted in Figure 7 and listed in Table 5. 
The IT profile shows excellent agreement with the results 
obtained from the joint Chandra- XMM fits. In particular, 
notice that it is the temperature values obtained for shells 
2C, 4C, and 6C that agree well with those obtained from 
the larger shells. The lower temperature values obtained in 
1C, 3C, and 5C indicate that the emission- weighted tem- 
peratures obtained in shells 1-3 are heavily weighted to- 
ward the outer edge of the shell. 

The temperatures for the 2T model in the smaller Chan- 
dra apertures also agree well with those obtained from the 
larger apertures. The value of Th in shell 1C is not shown 
since it has a large best-fitting value and is effectively un- 
constrained. Note also that results for a small number of 
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Table 4 

X 2 /dof for Selected Models Fitted Only to the Chandra Data 





Rin 


Rout 








Shell 


(arcmin) 


(arcmin) 


IT 


2T 


GDEM 


1C 


0.00 


0.25 


100.7/65 


94.1/63 


98.3/64 


2C 


0.25 


0.50 


130.9/85 


95.7/83 


118.2/84 


3C 


0.50 


1.00 


177.1/116 


122.4/114 


149.1/115 


4C 


1.00 


1.50 


196.5/120 


129.4/118 


141.5/119 


5C 


1.50 


2.00 


235.1/118 


159.1/116 


164.7/117 


6C 


2.00 


2.50 


170.1/116 


138.9/114 


148.0/115 



Note. — All models are 3D. These models also account for emission projected from shells outside shell 6C using the results obtained from 
fitting the XMM data outside shell 6C. 



Table 5 

Parameters for Selected Models Fitted Only to the Chandra Data 



Shell 


IT 

T 
(keV) 


2T 

T 
(keV) 


T h 
(keV) 


GDEM 
To a T 
(keV) (keV) 


1C 
2C 
3C 


0.61 ±0.04 
0.73 ±0.03 
0.75 ±0.01 


0.61 ±0.03 
0.70 ±0.03 
0.75 ±0.01 


2.1 ±0.6 

2.2 ± 1.2 


0.63 ± 0.11 
0.77 ±0.05 
0.76 ± 0.01 


0.15 ±0.14 
0.23 ± 0.08 
0.11 ±0.05 


4C 


0.79 ±0.01 


0.78 ±0.02 


1.58 ±0.35 


0.84 ±0.01 


U +0 ' 29 


5C 


0.88 ±0.01 


0.80 ±0.03 


1.45 ±0.13 


0.91 ± 0.03 


19+ 024 
u - ia -0.09 


6C 


0.91 ±0.01 


0.84 ±0.02 


1.27 ±0.23 


0.93 ± 0.03 


o io+ - 30 

u - ±u -0.03 



Note. — Models correspond to those in Table 4. Quoted errors using the symbol "±" are lc. Since the error ranges of try for shells 4C-6C 
are very asymmetric we quote the extreme values obtained from the 20 Monte Carlo error runs; i.e., these are effectively 95% confidence limits. 



5 10 15 20 5 10 15 20 




r (arcmin) r (arcmin) 



Fig. 7. — 3D radial temperature profiles (units - bottom: arcminutes, top: kpc) and la errors for (Left panel) IT and (Right panel) 2T 
models. The (blue) circles and solid diamonds refer to models fitted only to the Chandra data in shells 1C-6C. The (red) boxes and dashed 
diamonds refer to models fitted jointly to the XMM and Chandra data in shells 1-3. 



error simulations where the normalization of Th in shells 
2C and 3C is very small are not included in the error bars 
indicated in the figure. 

ft is clear that the 2T models fitted to the Chandra 
data alone indicate a range of temperatures at a given 



radius consistent with the 2T models fitted within wider 
shells. In shells 4C-6C the temperature is consistent with 
a constant value, Th between 1.3-1.5 keV. The tempera- 
ture difference within each shell is significantly larger than 
indicated by the radially varying single-phase model. For 
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example, in shell 5C where the ratio the emission measures 
for the cool and hot components is 2.3 (best fit), we have 
7h — T c = 0.65 ± 0.13 keV which is substantially larger 
(5a) than the ~ 0.03 keV spread in that shell expected 
from the single-phase hypothesis. 

Perhaps a clearer, more quantitative, indication of the 
temperature range within the shells is indicated by the 
GDEM model. The x 2 values and temperature param- 
eters are listed in Tables 4 and 5 are for the GDEM 
model. Similar to the results obtained from the joint 
XMM- Chan drafits in larger regions (§4.2.3), the x 2 values 
obtained for the GDEM model are generally intermediate 
between the IT and 2T models, though are clearly better 
than IT in all shells except 1C. The parameter values for 
the GDEM model obtained from the Chandra data alone 
are consistent with those obtained in the wider apertures 
(cf. Table 3) though with larger error bars. 

As we did in §4.2.3 we can use the GDEM model to 
test the hypothesis that the temperature profile of the 
hot gas represents a radially varying single-phase medium. 
Using the IT results for shells 3C-6C we would expect 
the following range of temperatures within shells 4C-6C: 
w 0.03 keV within 4C and 5C, and w 0.02 keV within 6C. 
Since the 95% confidence lower limit on ot is « 0.10 in 
these bins (Table 5), the lor temperature range implied 
by the GDEM model within each shell 4C-6C is at least 
0.20 keV - very much larger than than the 0.02-0.03 keV 
ranges expected from the single-phase hypothesis. (Con- 
sistent evidence for multiphase gas is indicated in shells 
2C and 3C, but the significance is only at the 2-3 o level.) 

We conclude that the improvement in the fits provided 
by the GDEM and 2T models, as well as the large implied 
temperature widths obtained for the Chandra data in the 
thinner shells, provides important additional evidence that 
a single-phase description of the hot gas in NGC 5044 is 
inadequate in the central regions. 

5. NON-RADIAL ANALYSIS 

5.1. XMM 

We performed a two-dimensional spectral analysis to de- 
termine whether the multiple components inferred from 
the analysis assuming spherical symmetry arise from az- 
imuthal temperature fluctuations within each annulus. We 
searched for azimuthal variations in the temperature and 
metal abundances using the following simple procedure; 
results obtained for the abundances are discussed in Paper 
2. Within the central CCD of the MOS images we defined 
a 5x5 array of equally spaced circular extraction regions of 
1' radius. Just outside of the central CCD, we defined 12 
equally spaced circular regions of 2' radius that surrounded 
the central CCD. Similar to the azimuthally averaged anal- 
ysis, to each region we fitted IT and 2T models modified 
by foreground Galactic absorption. Each model was fitted 
simultaneously to the MOS1 and MOS2 data projected on 
the sky; i.e., no deprojection. 

Overall, we find results consistent with the spherically 
symmetric analysis. The temperatures obtained from the 
IT model generally vary significantly only with distance 
from the center of the image. However, we notice a small 
asymmetry in the temperature distribution at a radius be- 
tween 2' - 3'. At this radius we find T w 0.9 keV for 
9 = 90° - 180° while T w 1.1 keV elsewhere at this ra- 



dius (6 is measured N-E.). This small azimuthal variation 
cannot account for the wide temperature distributions in- 
ferred from the 2T and continuous DEM models. 

In fact, when 2T models are fitted to the data, results in 
very good agreement with those obtained from the spher- 
ical analysis are obtained. The fits are improved over the 
IT models and values for T c and Th are obtained consistent 
with those presented for the 2T (2D) models in §4.2.2. 

5.2. Chandra 

We have performed a similar but independent two- 
dimensional analysis on the Chandra data. The higher 
spatial resolution of Chandra allows a more definitive test 
of the existence of large-scale multiple temperature gas 
components, by allowing regions as small as a few arcsec- 
onds in radius to be analyzed in principle. To directly 
compare with the azimuthally-averaged analysis, we have 
chosen circular regions which fall within annuli 2 and 3 
defined in Table 1. In Fig. 8 we present several represen- 
tative regions from this analysis overlaid on a smoothed 
Chandra image (we note that we have analyzed similar 
regions across the entire image, but do not find further 
variation beyond what we report below) . Six of these rep- 
resentative regions are 30" in radius, thus reaching the 
edge of the annulus they fall within. In addition, we have 
chosen three smaller regions 15" in radius which lie within 
the larger regions, and would fall within the thinner annuli 
3C and 5C (defined in Table 4). Some of the regions were 
placed intentionally in the areas of bright or faint emission 
visible in the smoothed image, which might be expected to 
exhibit the strongest spectral variations. As above, in each 
region we fitted IT and 2T models modified by foreground 
Galactic absorption. 

In Table 6 we present the results of the IT and 2T fits 
for the selected regions shown in Fig. 8. We again find 
results entirely consistent with the spherically symmetric 
analysis. For the IT model, the temperatures from regions 
lying within annulus 2 are clearly grouped about the value 
obtained for the entire annulus (cf. Figure 3), not devi- 
ating enough to explain the multi-temperature results of 
the annular analysis. In annulus 3, we see the same asym- 
metry as observed in the XMM non-radial analysis, with 
T w 0.8 keV for 9 = 90° - 180° while T w 1.1 keV else- 
where at this radius. Again, this variation cannot account 
for the wide temperature distributions inferred from the 
2T and continuous DEM models. 

The 2T model fits provide an impressive confirmation 
of this result. For regions 1A, 2, 3A, 6A, and 6B, the 
2T models are significantly improved over the IT models, 
and obtain values for T c and Th consistent with those fit 
to the entire annulus (§4.2.2). The strongest constraints 
are within annulus 3, which is expected since the emission 
measures of the hotter and cooler components are most 
similar there (§4.2.2; see Figure 4). For the four remain- 
ing regions, there is insufficient S/N to constrain the value 
of a second temperature component. In these regions (IB, 
3B, 4, 5), a second temperature component of similar tem- 
perature and relative emission measure strength as found 
in other regions can be added without significant effect on 
the fit. We emphasize that regions IB and 3B are simply 
subsets of their surrounding regions 1A and 3A, wherein 
a 2T model improved the fit. Presumably, in a longer ex- 
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Table 6 

Parameters from the Non-Radial Analysis of the Chandra, Data 







IT 




2T 










T 








tlc/llh 


Region 


X 2 /dof 


(keV) 


X 2 /dof 


(keV) 


(keV) 


(ratio) 


Annulus 2 














1A 


115.6/79 


0.82 ± 0.005 


81.7/77 


0.80 ± 0.009 


1.84 ±0.22 


3.0 ±0.6 


IB 


58.3/46 


0.76 ± 0.009 


53.8/44 


0.76 ±0.022 






2 


114.3/83 


0.78 ±0.005 


71.8/81 


0.74 ±0.015 


1.42 ±0.13 


3.0 ±0.7 


3A 


118.0/82 


0.81 ± 0.005 


79.6/80 


0.78 ±0.010 


1.52 ±0.19 


3.3 ±0.8 


3B 


42.5/42 


0.82 ± 0.010 


35.0/40 


0.82 ±0.015 






Annulus 3 










4 


47.6/41 


1.07 ±0.018 


45.9/39 


1.06 ±0.059 






5 


65.3/40 


1.12 ±0.031 


55.8/38 


1.04 ±0.067 






6A 


80.4/61 


0.77 ±0.013 


54.8/59 


0.73 ±0.017 


1.56 ±0.21 


2.9 ±0.8 


6B 


33.1/31 


0.75 ± 0.017 


26.6/29 


0.67 ±0.038 


1.14 ±0.19 


2.1 ± 1.1 



Note. — Models correspond to those in Table 4. Regions correspond to Figure 8. Regions with no entry for did not significantly constrain 
that parameter. Quoted errors using the symbol "±" are la. The final column is the ratio of emission measures of the cold and hot temperature 
components. 




Fig. 8. — Portion of the archival Chandra ACIS-S image located at the center of NGC 5044 adaptively smoothed using the CIAO task 
CSMOOTH. Heavy circles mark the limits of annuli 2 and 3. Overlaid are several regions used in our non-radial analysis (see §5.2). Region 
labels match Table 6. 



posure of the smaller regions, the same 2T result would 
obtain; i.e., while present, the 2T model is simply uncon- 
strained in the low S/N regime. 

Spatial variations in temperature could also be investi- 
gated with a hardness ratio map, an analysis which has 
been performed by Tamura et al. (2003). They do not 
find any significant azimuthal hardness variations, consis- 
tent with our non-radial analysis here. 

6. SYSTEMATIC ERRORS 

This section contains a detailed investigation of system- 
atic errors on the temperature measurements. Those read- 
ers who are not interested in these technical details can 
safely skip ahead to §7. 

6.1. Calibration 



6.1.1. XMM-EPIC vs. Chandra- ACIS 

We have examined possible systematic errors in the mea- 
surements of the temperatures arising from calibration dif- 
ferences between the XMM and Chandra CCDs. In Table 
7 we list the temperatures obtained from IT (2D) and 2T 
(2D) models fitted separately to the XMM and Chandra 
data; i.e., the MOS and pn data were fitted simultaneously 
while the Chandra data were fitted alone. We focus on 2D 
models so that the fits for a specific annulus are indepen- 
dent of results obtained from fits to adjacent regions at 
larger radii. 

Using the IT model the XMM and Chandra data in 
annuli 1 and 3 give values of Zp c that agree within 3% 
and are consistent within their 3er statistical errors. The 
temperatures of the 2T model are more uncertain, and 
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Table 7 

Comparison of Temperatures from XMM and Chandra 



IT 2T 







T 
















Ammlus 


XMM 


Chandra 


% 


XMM 


Chandra 


% 


XMM 


Chandra 


% 


1 


0.764 ±0.005 


0.738 ± 0.009 


3± 1 


0.747 ±0.006 


0.724 ±0.010 


3±2 


1.62 ±0.11 


2.08 ±0.31 


-28 ± 20 


2 


0.814 ±0.004 


0.794 ± 0.004 


3± 1 


0.790 ± 0.002 


0.782 ± 0.005 


1 ± 1 


1.46 ±0.04 


1.54 ±0.07 


-5 ±6 


3 


0.990 ± 0.004 


0.975 ± 0.005 


2± 1 


0.830 ±0.005 


0.835 ± 0.008 


-1±1 


1.36 ±0.03 


1.41 ±0.06 


-4 ±5 



Note. — All models are 2D to allow an independent comparison of the data sets in the central regions. The temperature is expressed in 
keV. "%" is the percent difference between the XMM and Chandra temperatures. 



the XMM and Chandra temperatures are consistent within 
their 1 — 1.5a errors. 

6.1.2. EPIC and ACIS CCDs vs. RGS Gratings 

We have compared results obtained from our analysis 
of NGC 5044 using the XMM and Chandra CCDs to 'the 
published results obtained for the XMM RGS obtained by 
Tamura ct al. (2003). First, by comparing the emission 
measures obtained within our extraction radii with the 
published emission measure for the RGS data by Tamura 
et al. (2003) and a separate analysis of the RGS data (T. 
Fang 2003, private communication) we conclude that the 
effective RGS aperture corresponds to a circular aperture 
with radius R ~ 30" — 40"; i.e., essentially our radius bin 
1. 

Tamura et al. (2003) find that a 2T (2D) model is a sig- 
nificant improvement over a IT model. In their analysis 
they (1) restrict their analysis to 0.44-1.55 keV, (2) only al- 
low the iron and oxygen abundances to vary, (3) use the so- 
lar abundance table of Anders & Grevesse (1989), and (4) 
use a plasma code essentially corresponding to the MEKAL 
code (Kaastra & Mewe 1993; Liedahl et al. 1995) in xspec. 
They obtain the following parameters for the 2T model: 
T c = 0.7 keV (no error quoted), T h = 1.07 ± 0.03 keV, 
Z Fc = 0.55 ± 0.05Z Q , and Z /Z Fc = 0.51 ± 0.06 in so- 
lar units. If we follow the same procedures (l)-(4) for the 
EPIC and ACIS data within radial bin 1 we obtain the 
following best-fitting values for the 2T model: T c = 0.68, 
21 = 1.08 keV, Z Fc = O.58Z , and Z /Z Fc = 0.46 in so- 
lar units. These results are in excellent agreement with 
Tamura et al.'s results for the RGS. 

We conclude that the RGS and the EPIC and ACIS 
CCDs give fully consistent results when the same models 
are fitted over the same energy ranges for each instru- 
ment. In particular, the RGS corroborates the improve- 
ment of a 2T (2D) model over a IT (2D) model. However, 
since they are unable to perform a deprojection analysis, 
Tamura et al. (2003) could not determine whether the ex- 
tra temperature component arises only from projection of 
gas from larger radii. (The consistency of the metal abun- 
dances is discussed further in Paper 2.) 

6.2. Plasma Codes 

We compared the results obtained using the APEC code 
to those obtained using the mekal code (Kaastra & Mewe 
1993; Liedahl et al. 1995) to assess the importance of dif- 
ferent implementations of the atomic physics and different 
emission line lists in the plasma codes. In every case ex- 
amined we found no qualitative differences between results 



obtained from each code; e.g., the fitted temperatures usu- 
ally agree to within ~ 5%. The x 2 values obtained for the 
multitcmperature models are also very similar. Some sig- 
nificant quantitative differences in \ 2 values are observed 
for the IT models, but there is no qualitative difference in 
the fits. For example, for shell 3 the IT (3D) mekal model 
gives x 2 — 1535.3 for 693 dof compared to x 2 — 1338.2 for 
the APEC code. Visual inspection of these fits and resid- 
uals reveals no noticeable differences between the two fits 
except for slightly more pronounced residuals in the Fe L 
region in the MEKAL fit. 

(We note that we explored the the validity of the as- 
sumption of ionization equilibrium using the vnei model 
in XSPEC. We found no improvement in the fit for a single- 
temperature model when allowing for departures from ion- 
ization equilibrium.) 

6.3. Bandwidth 

We explored the sensitivity of our results our default 
lower limit of the bandpass, E m i n = 0.5 keV. For com- 
parison we performed IT and 2T fits (both in 3D) with 
E min = 0.3 keV and E min = 0.7 keV. The fitted tem- 
peratures are generally consistent between models with 
-Emin between 0.3-0.7 keV. However, the x 2 values indi- 
cate that the improvement of the 2T model over the IT 
model decreases significantly as E m - m increases. In shell 
3, for example, we obtain x 2 /dof of 1427.9/774 (IT) and 
957.5/772 (2T) for E min = 0.3 keV , 1189.0/693 (IT) and 
836.0/691 (2T) for E min = 0.5 keV, 916.4/614 (IT) and 
730.5/612 (2T) for E min = 0.7 keV (each model is 3D). 
Similar behavior is observed for shells 1,3, and 4. 

This decreasing of the need for multitcmperature mod- 
els with increasing E m i n demonstrates that the residuals 
in the Fe L lines near 1 keV for IT models are not the 
sole driving force for the 2T and other multitcmperature 
models. This is significant since it implies that re- 
maining inaccuracies in the the Fe L lines in the 
plasma codes cannot be solely responsible for the 
improvement of the multitemperature models over 
IT models. Conversely, large bandwidth is seen to be 
imperative in the search for, and constraint of, multitem- 
perature models of the hot gas. 

6.4. Variable N H 

Since the lever-arm provided by the bandwidth below 
~ 0.7 keV is an important constraint on the multitempera- 
ture models, it should be expected that the fits may also be 
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Table 8 

IT and 2T Models with Variable Absorption 







"1 T 
1 I 




2T 






AJV H 




AN H 


Shell 


X 2 /dof 


(10 20 cm" 2 ) 


X 2 /dof 


(10 20 cm- 2 ) 


1 


688.7/450 


5.2 ±0.8 


637.8/448 


1.9 ±0.9 


2 


1059.8/688 


7.0 ±0.4 


916.7/686 


3.4 ±0.6 


3 


1120.9/692 


4.2 ±0.4 


828.2/690 


1.9 ±0.6 


4 


616.3/481 


4.0 ±0.7 


533.6/479 


1.3 ±0.8 


5 


510.3/453 


0.8 ±0.7 


485.2/451 


0.4 ±0.8 


6 


551.4/468 


0.8 ±0.8 


513.8/466 


0.4 ±1.0 


7 


561.2/514 


-1.8 ±0.8 


544.7/512 


-2.4 ±0.9 


8 


674.3/560 


0.5 ±0.7 


657.9/558 


0.8 ±1.0 







Note. — AJVh is the difference between the fitted absorption column density and the assumed Galactic value (5 X 10 20 cm 2 ). Both emission 
models are 3D while the absorption model is a conventional foreground screen. See §6.4 for description of the models. 



improved - to some extent - by allowing for intrinsic (con- 
tinuous) photoelectric absorption from cold gas. In Tabic 
8 we present results for the IT (3D) and 2T (3D) mod- 
els where we have allowed Nn of the foreground absorber 
component to be a free parameter. (We have examined a 
suite of cold absorber models and they give results simi- 
lar to those for the simple foreground-screen cold absorber 
model.) In shells 1-2 the IT model is improved signifi- 
cantly with variable absorption. However, in shells 1-4 the 
2T model is still clearly preferred over the IT model when 
allowing for variable absorption in each case. In fact, vari- 
able ./Vh actually improves the 2T model very little; i.e., 
the 2T model with Galactic absorption (Table 1) is itself 
superior to the IT model with variable N^m shells 1-4. 

The superior fits provided by the 2T (and other multi- 
temperature models like the PLDEM) with Galactic ab- 
sorption, and the relatively insignificant improvement to 
the multitemperature models obtained when intrinsic ab- 
sorption is allowed for, indicates to us that there is lit- 
tle motivation to consider intrinsic absorption models. 
Moreover, the variable-iVumodels imply large amounts 
of absorbing material in shells 1-4; e.g., taking AiVn = 
2.5 x 10 20 cm -2 in shells 1-2 leads to an absorbing mass, 
M abs w TOHAA^ H 7r(14kpc) 2 w 10 9 M Q assuming solar 
abundances and that the absorber is uniformly distributed; 
this is a lower limit if the abundances are sub-solar and the 
absorber is non-uniform. This large absorbing mass is al- 
most as large as that of the hot gas. Yet this amount 
of cold material, whether it be from optical line emitting 
gas or dust (e.g., Goudfrooij et al. 1994) or from atomic 
or molecular gas (e.g., Bregman et al. 1992; O'Dea et al. 
1994), has never been seen in NGC 5044 or other "cooling 
flow" galaxies. (We note that the fitted values of Nn for 
both the IT and 2T models at large radius are consistent 
with the Galactic value within their 1 — 2cr errors.) 

Since the multitemperature models with TVh = A^ al 
provide better fits within the central w 30 kpc, there is 
no obvious sharp absorption feature in the spectrum, and 
there is no evidence from observations in other wavebands 
for the large quantities of cold absorbing material implied 
by the fitted values of 7Vh, we do not take seriously the 
results obtained from the intrinsic cold absorber models. 

6.5. Background 



Since NGC 5044 is sufficiently bright the fitted tempera- 
ture values are quite insensitive to errors in the background 
normalization. For example, even if we do not subtract the 
background in the bounding annulus (where background 
is most important) the temperature we obtain for a IT 
model differs only by 10% from the background-subtracted 
value. Within the inner shells the background effect is neg- 
ligible; e.g., in shell 2 if the background is not subtracted 
we obtain a IT temperature that differs by 0.2% from the 
background-subtracted value. 

7. SPECTRAL ANALYSIS WITHIN A LARGE APERTURE 

Thus far we have fitted models to the XMM and Chan- 
dra spectra of a particular radial bin independently of 
the spectra in other radial bins. (Although a deprojected 
model for a specific radial bin does account for the pro- 
jected emission from exterior shells, the spectral fitting 
is not performed simultaneously with other shells.) Here 
we consider models fitted to the total XMM and Chan- 
dra spectral data accumulated within radial bins 1-3; i.e., 
within a circular aperture of radius, R = 2.5' (24 kpc). 
We focus on this central region where there is both Chan- 
dra and XMM data and evidence for multiphase gas. Our 
objective is to show that the preference for the multiphase 
models indicated by the fits to the individual radial bins 1- 
3 is rendered even more significant when the data in these 
bins are analyzed simultaneously in a single aperture. 

First we summarize the results (Table 9) of fitting IT, 
2T, and PLDEM models (all 2D for this exercise) to the 
spectral data within R — 2.5'. As expected, the IT model 
is an even worse fit than within the individual radial bins 
1-3. The best-fitting IT model shown in Figure 9 displays 
more pronounced residuals near 1 keV (i.e., characteristic 
of the "Fe Bias" ) than observed for the fit only to radial bin 
2 shown in Figure 5. The Chandra ACIS-S3 data have the 
most pronounced x 2 residuals in Figure 9 near 1 keV, but 
we note that the fractional residuals between the model 
and data are very similar for both the EPIC and ACIS 
CCDs. Allowing for variable -/Vh does not improve the IT 
fits significantly, despite the fact that a large excess value 
for the column density is indicated. 

However, the addition of a second temperature compo- 
nent (i.e., 2T model) provides a vastly improved fit with 
small residuals near 1 keV similar to the residuals at other 
energies (Figure 9). The PLDEM model fits nearly as 



One Temperature 



Two Temperatures 
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channel energy (keV) channel energy (keV) 

Fig. 9. — EPIC MOS1, MOS2, pn and ACIS-S3 spectra accumulated within a circular aperture of radius, R = 2.5' (24 kpc), fitted with 
(Left panel) a single temperature (IT) model and (Right panel) a two-temperature (2T) model (no deprojection is performed). In each case 
the APEC plasma model is used and the solar abundances are taken from Grevesse & Sauval (1998) which use the new (smaller) photosphcric 
value for the iron abundance. The models are the same as discussed in §4.2.1 and §4.2.2 for the spatially resolved analysis. That is, Galactic 
absorption is assumed and the following metal abundances are free parameters: O, Ne, Mg, Si, S, Fe and all other abundances are tied to Fe 
in their solar ratios. For the 2T model the abundances of each temperature component are tied together in the fits. 



Table 9 

Selected Results for Models Fitted Within a Large Aperture 



Model 


X 2 /d°f 


AJV H 
(10 2() cm- 2 ) 


T c 
(keV) 


T h 
(keV) 


± mm 

a (keV) 


2"max 2"min 

(keV) 


IT, Galactic N H 
IT, variable JV H 
2T, Galactic N H 
2T, variable N H 
PLDEM, Galactic TVh 


4197.4/869 
3879.5/868 
1356.8/867 
1313.0/866 
1373.0/867 


4.5 ±0.3 
3.0 ±0.6 


0.897 ±0.002 
0.878 ± 0.002 
0.792 ±0.002 
0.779 ±0.006 


1.41 ±0.02 
1.28 ±0.04 


-3.0 ±0.1 0.681 ±0.005 


1.5 ±0.1 



Note. — The (2D) models are fitted to the accumulated XMM EPIC MOS1, MOS2, pn and Chandra ACIS-S3 spectral data within a 
circular aperture of radius 2.5' (24 kpc). A/Vjj is the difference between the fitted absorption column density and the assumed Galactic 
value. The emission measures for the 2T models are, in the XSPEC units discussed in the caption to Figure 4: norm c = 6.7e-3 ± le-4 and 
norm h = 3.7e-3 ± le-4 for 2T, Galactic JV H and norm c = 8.1e-3 ± 3e-4 and norm h = 4.6e-3 ± 4e-4 for 2T, variable A^h- 



Table 10 

Fits to Simulated Data of Radially Varying Models Within a Large Aperture 



Model 


Sim #1 


X 2 /dof 
Sim #2 


Sim #3 


IT, Galactic N H 
IT, variable iV H 
2T, Galactic N H 


1809.8/839 
1805.3/838 
865.6/837 


2469.2/855 
1309.1/854 
1105.4/853 


3731.6/861 
3611.3/860 
980.0/859 







Note. — The models are fitted to simulated XMM EPIC MOS1, MOS2, pn and Chandra ACIS-S3 spectral data accumulated within a 
circular aperture of radius 2.5' (24 kpc). "Sim #1-3" refer to simulations of models obtained from the spatially resolved analysis in radial bins 
1-3 following our discussion in §7. 



well as the 2T model, and the quality of both of these 
fits is about as good as could be hoped for considering 
the simplicity of these models and that we know the spec- 
tral properties do vary with radius within the R = 2.5' 
aperture. (Allowing for variable A% also provides negli- 
gible improvement for the multitemperature models even 
though a large fitted value of excess TVh is indicated - the 
2T result is shown in Table 9.) 

Now we wish to compare these results to those that 
should have been obtained if the real data were actually 
described by the radially varying IT and 2T models in 



shells 1-3 obtained in §4.2.1 and §4.2.2. To perform this 
comparison we follow the general approach discussed in 
§5.3 of Buote (1999). We take the IT (3D) models in 
shells 1-3 obtained from the real EPIC and ACIS data 
(§4.2.1) and simulate EPIC and ACIS spectra appropri- 
ate for the NGC 5044 observations using the fakeit rou- 
tine in xspec. These simulated spectra also contain the 
projected emission from external shells obtained from our 
deprojection analysis. The resulting simulated source and 
background pha files for radial bins 1-3 arc then summed 
(separately for each detector) and analyzed in the same 
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way as done in Table 9. We refer to this simulated radially 
varying IT model within R = 2.5' as "Sim #1" in Table 
10. "Sim #2" and "Sim #3" in Table 10 are prepared in 
the same way as "Sim #1" except they refer respectively 
to the IT (3D) model with variable N H (§6.4) and 2T (3D) 
model with Galactic A H (§4.2.2). 

Fits to Sim #1 using IT models give x 2 ~ 1800 which 
are much smaller than the values of \ 2 ~ 4000 obtained for 
the real data in Table 9. The 2T model is also a somewhat 
better fit than for the real data and provides a formally ac- 
ceptable fit. Although similar to the real data the 2T fit is 
superior to the IT fit, and there is negligible improvement 
when allowing for variable absorption, the much smaller 
values of x 2 f° r the IT models fitted to Sim #1 rule out 
the radially varying IT models with Galactic absorption. 

Qualitatively different results are obtained for fits to 
Sim #2. Similar to Sim #1 the fitting the IT model with 
Galactic Anhas yields a much lower value of \ 2 than ob- 
tained for the real data. However, the IT model with 
variable Nn provides a very large improvement in the fit 
- almost as good as the 2T model with Galactic Ah- This 
behavior is totally inconsistent with the results obtained 
from the real data in Table 9, and we conclude that the 
radially varying IT model with variable Ah is not a viable 
description of the spectral data within R — 2.5'. 

In contrast, the \ 2 values obtained for fits to Sim #3 
are very similar to those obtained for the real data in Ta- 
ble 9. Since only fits to Sim ^3 can reproduce the results 
of fitting the real data, this is strong evidence that the 
multitemperature models (in this case the radially vary- 
ing 2T model, but this also applies to the PLDEM model) 
are required within the central R = 2.5' (24 kpc) of NGC 
5044. 

8. DISCUSSION: TWO-PHASE MODEL 

Our spectral deprojection analysis of the XMM and 
Chandra data of NGC 5044 indicates that a single-phase 
description of the hot gas is inadequate for r < 30 kpc. 
Of the simple multiphase models we considered, the 2T 
and PLDEM models provide the best fits in the central 
regions. Each of these models describes a limited multi- 
phase plasma where the cooler temperature components 
(<~ 0.7 keV) dominate for r < 10 kpc, the contributions 
of cooler and hotter (~ 1.4 keV) components are simi- 
lar for r ~ 20 — 30 kpc, while at larger radii the hotter 
components dominate so that the gas is consistent with a 
single-phase medium. 

These results for the 2T model are very similar to those 
obtained from an analysis of the XMM data of the group 
NGC 1399 (Buote 2002). The temperature of the ex- 
tended hotter component, Th <~ 1.4 keV, is consistent 
with the virial temperature of a surrounding group of mass 
ss 10 13 M Q , whereas the temperature of the centrally con- 
centrated cooler component, T c <~ 0.7 keV, is similar to the 
kinetic temperature of the stars. These parameters sug- 
gest a physical association of the hotter component with 
the ambient group gas and the cooler component with stel- 
lar ejecta from the dominant central galaxy. 

A potential problem with this scenario is that it might 
be expected that these phases should mix very rapidly for 
a system in equilibrium. However, the sharp edge in the 
surface brightness and isophote center offset for R 5' — 6' 



discussed in §3 are suggestive of a "cold front" such as 
has been observed in several clusters with Chandra (e.g. 
Markcvitch et al. 2002). The surface brightness edge in 
the MOS image of NGC 5044 near R w 50 kpc to the NW 
essentially divides the regions where the cooler and hotter 
gas phases are most prominent. 

Although the cold front can explain the radial transi- 
tion from cooler to hotter dominance, the phases clearly 
co-exist over a large range in radius. For this to occur 
we expect the phases to be in pressure equilibrium. We 
have calculated the volume filling factor (/ c ) of the cooler 
component of the 2T model required to maintain pressure 
equilibrium between the two phases. The result is plotted 
in Figure 10. 

For r > 30 kpc we have f c <~ 0.1 indicating that the 
cooler gas occupies only a small fraction of the shell vol- 
ume. In contrast, for shells 2-3 (r <~ 5 — 25 kpc) we have 
f c w 0.5 meaning that each phase occupies half the vol- 
ume. In shell 1 (r < 5 kpc) the value shown (/ c 0.3) is 
actually a lower limit because the value of Th is overesti- 
mated due to contamination from discrete sources (§4.2.2). 

9. CONCLUSIONS 

The spectral deprojection analysis of XMM and Chan- 
dra data favors a two-phase (2T) or limited multi-phase 
medium (PLDEM) within the central r ~ 30 kpc of NGC 
5044. The cooler component in the 2T models has a tem- 
perature T c ~ 0.7 keV similar to the kinetic temperature 
of the stars in the central galaxy NGC 5044, and the hotter 
component has a temperature Th ~ 1.4 keV characteris- 
tic of the massive ~ 10 13 M© dark halo of the surround- 
ing galaxy group. Nevertheless, both temperature com- 
ponents appear at every radius < 30 kpc. In spite of the 
similarity of the hot phase temperature Th and the group 
virial temperature at all radii, it is likely that gas at small 
radii with temperature Th is heated by a central AGN. 
Some additional heating at large radii could arise from 
the energy associated with establishing the cold front near 
r ~ 30 kpc, as suggested by the sharp edge visible in the 
EPIC MOS and pn images. 

The need for two discrete temperatures in the 2T models 
cannot be attributed solely to the Fe-L lines at <~ 1 keV, 
which still may be uncertain in the plasma codes, but is 
also required by the spectrum below 0.7 keV. Our com- 
parison of results using the APEC and MEKAL codes in §6.2 
indicates that even large differences in the plasma codes 
have a small effect at the moderate resolution of the EPIC 
and ACIS CCDs. 

As an alternative, and possibly more physically plausi- 
ble model, we have shown that a continuous, but limited, 
range of temperatures in each spherical shell (PLDEM) 
can explain the NGC 5044 data for r < 30 kpc as well 
as the 2T model. However, the range of temperatures 
required in each spherical shell exceeds the radial tem- 
perature variation of the best-fitting single-phase models 
across the same shell. For either type of thermal model 
there is no evidence at any radius for gas at temperatures 
< 0.7 keV. These results are very similar to our previous 
2T models of NGC 5044 using ASCA data (Buote 1999). 
However, within r <~ 30 kpc the ASCA data lacked the 
spatial resolution of XMM and Chandra and could not dis- 
tinguish between a single-phase medium in which the gas 
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Fig. 10. — Volume filling factor of the cooler component required to maintain pressure equilibrium between the cooler and hotter phases 
in the 2T (3D) model obtained from simultaneous fitting of the XMM and Chandra data. Due to the expected contribution of emission from 
discrete sources, the value in the central bin is quite uncertain and is consistent with unity. The units are arcminutes on the bottom axis and 
kpc on the top axis. 



temperature varies with radius, a two-phase (2T) medium, 
or multi-phase gas at every radius. 

The remarkable irregularities visible for r < 10 kpc in 
the Chandra image of NGC 5044 (Figure 1) support the 
notion of fluctuations in the gas density and, by implica- 
tion, also in the gas temperature. In pressure equilibrium 
the ratio of X-ray emissivities in the two phases should be 
approximately, [{n c ) c /(n c ) h ] 2 oc (T h /T c ) 2 w 5, which may 
be sufficient to account for the conspicuous surface bright- 
ness fluctuations in Figure 1. Within the central annulus 
(R = 0.5' sa 5 kpc), which includes nearly half the opti- 
cal image of the NGC 5044 galaxy, the X-ray spectrum is 
further complicated by hard ~ 10 keV bremsstrahlung ra- 
diation from X-ray binary stars, so the precise range of gas 
temperatures is less accurately determined in this region. 

The 2T or limited multi-phase thermal properties of 
NGC 5044 are very similar to those of the bright galaxy 
group NGC 1399 revealed by XMM (Buote 2002). A 2T 
model is also preferred by Molendi (2002) in his XMM 
analysis of M87. However, Molendi argues that the tem- 
perature range Th — T c implied for each annulus by his 2T 
model for M87 is consistent with a radially varying single- 



phase temperature except in regions where the X-ray im- 
age is clearly distorted by interaction with the central ra- 
dio source. Nevertheless, the pattern that is emerging from 
observations of NGC 5044 and NGC 1399 may differ sub- 
stantially from the gasdynamical models of groups and 
clusters constructed by Brighenti & Mathews (1999, 2002). 
In the gasdynamical models relatively cool ejecta from stel- 
lar mass loss mix thermally on small scales with hotter am- 
bient group or cluster gas to reproduce the single-phase ra- 
dial temperature profiles typically observed in groups and 
clusters. Instead, our observations suggest a limited range 
of thermal phases that are incompletely mixed at every ra- 
dius, which, if viewed as a single-phase gas, can reproduce 
the same typical "cooling flow" thermal profile but with a 
reduction in the quality of the spectral fit. 
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results of his preliminary analysis of the RGS data of 
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